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Abstract 

We compute, within the Schwinger-boson scheme, the Gaussian-fluctuation 
corrections to the order-parameter stiffness of two frustrated quantum spin 
systems: the triangular-lattice Heisenberg antiferromagnet and the Ji — J2 
model on the square lattice. For the triangular-lattice Heisenberg antiferro- 
magnet we found that the corrections weaken the stiffness, but the ground 
state of the system remains ordered in the classical 120° spiral pattern. In the 
case of the Ji — J2 model, with increasing frustration the stiffness is reduced 
until it vanishes, leaving a small window 0.53 ^ rj ^ 0.64 where the system 
has no long-range magnetic order. In addition, we discuss several method- 
ological questions related to the Schwinger-boson approach. In particular, we 
show that the consideration of flnite clusters which require twisted boundary 
conditions to flt the infinite-lattice magnetic order avoids the use of ad hoc 
factors to correct the Schwinger-boson predictions. 

I. INTRODUCTION 

Two-dimensional (2D) frustrated quantum spin systems are a fascinating subject that 
has been studied for more than 20 years now. Since the seminal work by Fazekas and Ander- 
son on the S = ^ triangular- lattice Heisenberg antiferromagnet (TLHA), |l| the existence 
of non-magnetic ground states in these systems has been strongly debated in the literature. 
In the last years, this problem was somehow linked to more general questions concerning 
spin-liquid states and their possible connections to high-Tc superconductivity. In this con- 
text, the square-lattice Heisenberg antiferromagnet with positive first- and second- neighbor 
interactions, the so-called J1—J2 model, has been much studied. Unlike the TLHA which 
is frustrated because of the lattice topology, in the Ji — J2 model frustration is due to the 
competing couplings. 

Regarding the TLHA, although in general there is a growing conviction that this model 
displays the characteristic 120° spiral order in its ground state (perhaps with an impor- 
tant reduction of the classical | value), using different methods some authors found a 
situation very close to a critical one or no magnetic long-range order at all. [^,0 We have pre- 
viously performed a study of the TLHA — including first- and second-neighbor interactions — 



using the rotational invariant Schwinger-boson approach in a mean-field (saddle-point) ap- 
proximation. In particular, for the ground-state energy this method showed a remarkable 
agreement with exact numerical results on finite lattices (at the time of this work only re- 
sults for the 12-site cluster were available). Furthermore, the local magnetization in the 
thermodynamc limit was predicted to be slightly larger than the (harmonic) spin-wave re- 
sult. More recently, new finite-size results and 1 /^-corrected spin- wave calculations 
PJTm became available, which lead us to reconsider this model. In this work we present the 
one-loop corrections to the saddle-point results obtained in , along the lines developed in 
rT|,T2|. We have computed the ground-state energy and the spin-stiffness tensor in order 



to assess the existence or not of magnetic order, and we have compared these results with 
exact values on finite lattices and spin-wave results. 

As for the J1—J2 model, in previous publications we have considered its behavior 

on finite clusters that do not frustrate the short-range Neel and collinear orders present in 
this system. In particular, we computed the stiffness p to determine the range of 77 = J2/J1 
where the magnetic order is destroyed by the combined action of quantum fluctuations and 
frustration. However, it turned out to be that the value of p obtained at saddle-point order 
was nearly a factor 2 smaller than exact results on small clusters. This was conflrmed by 
considering the large-S* limit of the Schwinger-boson predictions, which differed from the 
classical results by exactly this factor. The same happened, now for all values of 5, in the 
2-site problem, which can be worked out exactly. After these checks, in we simply 

added an ad hoc factor 2 to correct the saddle-point result for p. Although the need for this 



kind of factors has been stressed since the flrst works on Schwinger bosons, [|1^ their use is 
highly unsatisfactory. 

In the course of the present study on clusters of the triangular lattice, again the large-S" 
Schwinger-boson prediction for p did not reproduce the classical result. The situation was 
here even worse than for the square lattice, since now the rotational invariant Schwinger- 
boson approach could not produce concrete predictions for the stiffness parallel and per- 
pendicular to the plane where the order parameter spirals. In an attempt to resolve this 
last problem, we followed the idea in and considered clusters with appropiate twisted 

boundary conditions to flt the 120° spiral order in the a;?/-plane. In this case, since the 
boundary conditions break the rotational invariance explicitly, one is able to compute the 
parallel and perpendicular stiffness separately. Moreover, we found that now the large-S* 
predictions for these quantities have the correct behavior and no ad hoc factors are required. 
This unexpected result prompted us to reinvestigate the J1 — J2 model on clusters which re- 
quire the use of antiperiodic boundary conditions to avoid frustrating the Neel and collinear 
orders. Again in this case the results have the expected large-S" behavior, and there is no 
need for correction factors. These flndings point to a complex interplay between rotational 
invariance and the unphysical enlargement of Fock space in the Schwinger-boson approach. 
Fortunately, working on appropiate clusters the results for both the TLHA and the Ji — J2 
model behave consistently, as we show below. 



II. ORDER-PARAMETER STIFFNESS AND THE SCHWINGER-BOSON 

APPROACH 

Let's consider a general Heisenberg model. 
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H = lj:Jir-r')S^.Sr, (1) 

r,r' 

where r, r' indicate sites on a two-dimensional lattice. On finite clusters generated by 
the translation vectors (a = 1, 2) we impose arbitrary boundary conditions S^+t^ = 
7^-;j($Q,)S'r, where 7?.^($q,) is the matrix that rotates an angle $q, around some axis n (notice 
that we are using boldface (arrows) for vectors in real (spin) space). By performing local 
rotations Sr — ^ lZi;i{Or)Sr of angle 6'r = Q ■ r, the Hamiltonian (|l|) becomes 



if = - ^ J(r — r') |(S'r ■ n)(SV' ■ n) + cos6'rr' Sr ■ S^' — (5V ■ n)(SV' ■ n) + sin6'rr'(>S'r x SV') ■ nj 

r,r' 

(2) 

where 9rr' = Q ■ (r' — r). In this way, with the choice Q ■ = $o the boundary conditions 
become the standard periodic ones Sr+T^ = Sr- 
We define the (T = 0) stiffness tensor p;^ by 



^0 d'EcsiQ) 



(3) 

Q=0 



where Eqs is the ground-state energy per spin and = Q ■ e„ (a= 1, 2) are the twist angles 
along the basis vectors Bq,. Then, by using second-order perturbation theory it is simple to 
prove that pi^^ = T"^ -|- J"^, where 

Tf = i'l E [5^0 ■ Sr-{So ■ n){Sr -n)]) , (4) 



and 



Here = | Er ^(r)r"(S^o x SV) • n, P = 1 - |GS )( GS| , and r = r°e„. For further use, 
it is convenient to consider also the rotational average 'p^y of the stiffness tensor for n in 

the xy-p\ane. With n = (sin 6* cos 0, sin 5 sin 0, cos 6*) and 6 = using J^^ ^n^n^ = ^ 
= 1,2) one obtains p^y = lip^ + py)- Analogously, using J ^n^n^ = ^ = 1,3), 
for the corresponding average in all directions in spin space one obtains p = ^(pj + p^ + Px)- 
In the case of classical spins, the only contribution to the stiffness comes from the T-term 
(I). For quantum spins, besides the fact that the calculations are much more complicated, 
there is also a very important conceptual difference with the classical case. This is more 
clearly explained for the simple square-lattice Heisenberg model, where from the Lieb-Mattis 
theorem we know that the exact ground state must be an isotropic singlet 5 = (there is 
numerical evidence |P that this theorem — in principle valid for bipartite lattice — might be 
conveniently extended to the triangular and other frustrated lattices). Then, if we approach 
the infinite-lattice behavior using a technique which respects the rotational invariance — 
like, for instance, the Lanczos method on small clusters — , we only have access to the 



GS 
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stiffness average p. On the other hand, if we use approximate methods that assume an 
order-parameter symmetry breaking — hke ordinary spin-wave theory, for instance — , then 
both the parallel p\\ and perpendicular stiffness tensors can be computed. There is yet 
a third possibility, the use of the Schwinger-boson technique, which keeps the rotational 
invariance on finite lattices and develops a spontaneous symmetry-breaking (through a Bose 
condensation of the Schwinger bosons) in the thermodynamic limit. We will consider this 
last method since, in addition to providing us with a fairly reliable computational tool, there 
are several interesting methodological questions which can be discussed. 

First, we represent spin operators in terms of Schwinger bosons: \T2\ Si = |aj.(j.cij, where 
a; = (a||,a||) is a bosonic spinor, a is the vector of Pauli matrices, and there is a boson- 
number restriction J2a '^la'^io- = 25* on each site. Thus, we can write the partition function 
for the Hamiltonian as a functional integral over boson coherent states, which allows a 
saddle-point expansion to be performed. Following the calculations in [|T^], which we do not 
reproduce here, we obtain the fluctuation-corrected ground-state energy as the sum of the 
saddle-point and one-loop contributions, ii^Gs(Q) = -£'o(Q) + -E'i(Q)- We evaluate £'gs(Q) 
on finite lattices in order to avoid the infrared divergencies associated to Bose condensation, 
which signals the appearance of magnetic long-range order in the ground state. Numerical 
differentiation of -Egs(Q) according to gives finally the spin stiffness. 

In what follows we will use this method to study the order-parameter stiffness of both 
the TLHA and the J1 — J2 model. 



III. TRIANGULAR-LATTICE HEISENBERG ANTIFERROMAGNET 

For the triangular lattice, finite clusters with the spatial symmetries of the infinite lattice 
correspond to translation vectors Ti = {n + m)ei + me2, T2 = nei + {n + m)e2, where 
ei = (a, 0), e2 = (— |a, ^a) (a is the lattice spacing). The number of sites in these 
diamond-shaped clusters is given by = -|- -|- nm. When 2n + m and n — m are 
multiples of 3, which corresponds to clusters with N = 3p {p integer), periodic boundary 
conditions do not frustrate the 120° spiral order. Otherwise, N = 3p+l and one has to use 
twisted boundary conditions with angles |7r(2n + m) and |7r(r2 — m) along the Ti and T2 
directions, respectively. 

Let's consider first the (by far) easiest case of classical spins. Choosing the proper 
boundary conditions as explained above, the classical energy is minimized by a configuration 
Sr = S{cos9r, sin6'r, 0). Here Or = Qf^irai ■ with the magnetic wavevector Qf^irai = (|^! 0). 
Since we took the spins to lie on the xy-pla.ne, it is necessary to consider both the parallel 
stiffness p\\ = p^ under a twist around z, and the perpendicular stiffness p± = p^y for 
twists around n versors lying on this plane. As mentioned before, for classical spins the 
only contributions to py, p± and p come from the T term (^, and one always obtains the 
tensorial structure p"^ = |(1 + 6"f^)p with py = JS'^a?, p± = |p|| and p = |p||. These last 
two results are consecuences of the identities p^y = \{p^ + p^ and p = |(p2 + p^ + p^), plus 
the fact that p^ = |p|| and p^ = |p|j for the classical spiral arrangement we choosed. 

In order to correct these results by considering the quantum nature of the spins we have 
followed the procedure discussed in the previous section. On the N = 3p clusters with 
periodic boundary conditions our approach is rotationally invariant and we only have access 
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to p. Since now py and p± are unrelated, the theory cannot produce concrete predictions 
for them from p. On the other hand, we found that the large-S* Schwinger-boson prediction 
for this last quantity is exactly 4/3 smaller than the corresponding classical result. As 
mentioned in the introduction, in an attempt to resolve the former problem we followed the 
idea in P,p!0|, and considered the N = 3p + l clusters with the appropiate twisted boundary 
conditions to fit the 120° spiral order in the xy-plane. In this case, since the boundary 
conditions break the rotational invariance explicitly, one is able to compute the parallel and 
perpendicular stiffness separately. Furthermore, we found that now the large-S* predictions 
have the correct behavior and no ad hoc factors are required. 

In Fig. 1 we present the results for the energy of the N = 3p clusters, together with 
results from exact diagonalization studies As can be seen, for N = 12,36, after 

the inclusion of Gaussian fluctuations our theory gives predictions very close to the exact 
ones. However, for N = 21, 27 the agreement is not so good. This can be explained by 
noticing that for odd number of sites the true ground state has total spin S = 1/2, while the 
approximate Schwinger-boson wave function is always a singlet (rotational symmetries are 
broken only in the thermodynamic limit by the boson condensate). The difference between 
even and odd is also apparent in the exact results, since these do not line up according to 
the expected (spin- wave) scaling. However, following after substracting the top inertial 
effects in the exact results for the odd-A^ clusters, our results and these corrected values 
are in very good agreement (see Fig. 1). This problem becomes less important for larger 
clusters. The inset in Fig. 1 shows the thermodynamic-limit extrapolation on large lattices 
of both the N = 3p and N = 3p + 1 types (notice that for the N = 3p clusters the linear 
trend has a slope rather different than what one would obtain from the consideration of the 
smallest clusters). As discussed above, the A^ = 3p + 1 clusters require twisted boundary 
conditions to avoid frustrating the spiral order. The main consecuence of this is the absence 
of rotational symmetry in the ground state, since the boundary conditions choose a plane 
of magnetization (the xy-plane in this case). Consequently, the scaling of the energy to the 
thermodynamc limit has a quite different slope than the corresponding to A^ = 3p clusters, 
although, of course, the final value for the ground-state energy per site Eqs — —0.5533 is 
the same in both cases. This happens, as expected, with the results at saddle-point order 
and also after the inclusion of Gaussian fluctuations. 

The finite-size scaling of the parallel stiffness is shown in Fig. 2. In this figure we plot 
both the saddle-point and fluctuation-corrected results, which run almost parallel to each 
other. The first-order spin- wave result of is also given for comparison. Notice that, 
in spite of the different slope, the extrapolated values do not differ much. As a further 
comparison, in Fig. 3 we plot the renormalization factor between quantum and classical 
results for the parallel stiffness on small lattices, as given by different techniques. As can 
be seen from this figure, the introduction of Gaussian-fluctuation corrections to the saddle- 
point Schwinger-boson results gives values closer to the exact ones than those of first-order 
spin-wave theory. 



IV. J1-J2 MODEL 

For the square lattice, finite clusters of A^ = ^/n'^ + m? x i/n^ + m? sites correspond to 
translation vectors Ti = ne.i + me2, T2 = — mei + ne2i where ei = (a,0), 62 = (0,a). 
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If n, m are even, the periodic boundary conditions do not frustrate the short-range Neel 
(QNeei = (~) ~)) colhnear (Q°oii = (-,0)) orders present for 7]= J2/J1 -C 1 and 77 ~ 1 
respectively. For n and/or m odd one has to impose antiperiodic boundary conditions to fit 
both magnetic patterns to the cluster shape. 

Let's assume that the classical energy is minimized by a spin configuration = 
^(cos^^r, sin6'r, 0). For simplicity we consider first the square lattice with nearest-neighbor 
interactions only, so that 9r = QnccI ' ^- Then, from (|p one obtains the tensorial structure 
p"^ = pd"^, where, like for the triangular lattice, py = JS'^a?, p± = and p = 
Notice however that for colhnear spin arrangements like the Neel order there is properly 
no perpendicular stiffness. In fact, the value of p± is obtained from p^y = i(pj + p^) and 
/>2 = 0, = Pll for vectors pointing in the ±x directions. In the same way, the value of p is 
a consequence of the identity p = |(pj + + Pj)- 

In previous publications ||Tl| , p!2| we have considered finite clusters with periodic boundary 
conditions. Then, we determined p and, from this quantity, py = |p. However, it turned out 
to be that the results obtained at saddle-point order for py were nearly a factor 2 smaller (or, 
like in the triangular lattice, the results for p nearly 4/3 smaller) than exact results on small 
clusters. As pointed out in the introduction, this was confirmed by considering the large-S* 
limit of the Schwinger-boson predictions, which differed from the classical results by exactly 
these factors. Moreover, the same happened for all values of S in the 2-site problem, which 
can be worked out exactly. After these checks, in |lTl| , [T2| we simply added an ad hoc factor 
2 to correct the saddle-point result for py. The need for this kind of factors is related to a 
subtle interplay between rotational invariance and the relaxation of the local boson-number 
restriction, a claim that can be substantiated by the following argument. At saddle-point 
order the energy Egs(Q) = | Er </(r)[5r (Q) - ^r(Q)], where 5r(Q) and Ar{Q) are the 
order parameters for the spin-spin interaction in the ferromagnetic and antiferromagnetic 
channels respectively. At the same order in the calculations, the exact (local) constraints 
would force the identity B^{Q) + Al{Q) = S*^, which shows that both order parameters 
should contribute in the same way to the stiffness (^. Moreover, the same identity indicates 
that + = -{Ar gqgQ +Br q^q^ ) . However, when the constraints are imposed 

on average the order parameters behave independently, the identity is violated, and this 
last relationship is not fulfilled (the l.h.s. vanishes). On the other hand, on clusters with 
antiperiodic boundary conditions, which break the rotational invariance, this relationship 
between derivatives does come right and leads to the correct value for the stiffness, as shown 
below. 

In Fig. 4 we plot the spin stiffness of the J1—J2 model on large clusters with antiperiodic 
boundary conditions (we considered rj = and 0.5 as examples). As expected, the points line 
up according to the scaling behavior py ~ N~^^'^. This extrapolation to the thermodynamic 
limit for several values of t] produces the result shown in Fig. 5. There we see a small window 
(~ [0.53,0.60]) where the stiffness vanishes and there is no magnetic order (compare with 
the bottom panel in Fig. 3 of [0]). The inset shows the scaling behavior of the point 
'7Neei where the Neel-order stiffness vanishes. Actually, the region without long-range order 
extends at least up to 77 ^ 0.64, since, as shown in [l^, for 0.60 ^ 77 ^ 0.64 the short-range 
antiferromagnetic order has lower energy than the long-range colhnear one. 



Let's consider also the scaling of the ground-state energy with A^, Eq^^Q^ 



Neel) 



-E'gs(Qncci) ~ a/N^^"^, along the series of clusters with periodic and antiperiodic bound- 
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ary conditions. The inset in Fig. 6 shows that both series extrapolate to the same value 
(E'gsIQncsci) — —0.5015 for rj = 0.5, within numerical errors), although with quite different 
slopes a. It is of some interest to discuss the behavior of a, which, for periodic boundary 
conditions, is proportional to the spin- wave velocity (or its anisotropy average) Work- 
ing on clusters with periodic boundary conditions, in our previous publication |12[ we found 
an anomalous behavior of a near the point where the collinear order is melted (see Fig. 3 
in [0). In the present work, for clusters with antiperiodic boundary conditions, a has a 
behavior more in line with a priori expectations (see Fig. 6). 



V. CONCLUSIONS 

We have considered the Gaussian-fluctuation corrections to the spin stiffness of two 
frustrated quantum spin systems of interest: the TLHA and the Ji — J2 model. For the 
TLHA we found that the fluctuation corrections weaken the stiffness, but the ground state 
of the system remains ordered in the classical 120° spiral state. In the case of the Ji — J2 
model, with increasing frustration the stiffness is reduced until it vanishes, leaving a small 
window 0.53 ^ rj 0.64 where the system has no long-range magnetic order. 

In the course of this investigation we discussed several methodological questions related 
to the Schwinger-boson approach we used. In particular, we showed that the consideration 
of clusters which require twisted boundary conditions to fit the magnetic orders avoids the 
use of ad hoc factors to correct the Schwinger-boson predictions. This fact points to a 
subtle interplay between rotational invariance and the relaxation of local constraints in this 
approach. 

Finally, it is interesting to notice that for both the square and triangular lattices with only 
nearest-neighbor interactions the fluctuation-corrected spin stiffness scale almost parallel to 
the saddle-point results (see Figs. 2 and 4), although the corrections have different sign. 
On the contrary, for large frustration in the square lattice the one-loop result have a very 
different slope than the zero-order one, which ultimately leads to disordering the system. 
This difference reflects the distinct nature of the order-parameter fluctuations in this highly 
frustrated system. 



7 



REFERENCES 



[1] P. Fazekas and P. W. Anderson, Phil. Mag. 30, 423 (1974). 

[2] A. V. Dotsenko and O. P. Sushkov, Phys. Rev. B 50, 13821 (1994) and references 
therein. 

[3] B. Bcrnu, C. Lhuilher, and L. Pierre, Phys. Rev. Lett. 69, 2590 (1992); 

[4] D. A. Huse and V. Elser. Phys. Rev. Lett. 60, 2531 (1988); R. Deutscher and H. U. 

Everts, Z. Phys. B 93, 77 (1993). Phys. Rev. Lett. 60, 2531 (1988). 
[5] B. Bernu, P. Lecheminant, C. Lhuilher, and L. Pierre, Phys. Rev. B 50, 10048 (1994). 
[6] H. Nishimori and H. Nakanishi, J. Phys. Soc. Jpn 57, 626 (1988); 58, 2607 (1989); 58 
3433 (1989); N. Eltsncr, R. R. P. Singh, and A. P. Young, Phys. Rev. Lett. 71 , 1629 
(1993); K. Yang, L. K. Warman , and S. M. Girvin, Phys. Rev. Lett. 70, 2641 (1993). 
[7] P. W. Leung and K. Runge, Phys. Rev. B 47, 5861 (1993). 
[8] C. J. Gazza and H. A. Geccatto, J. Gondens. Matt. 5, L135 (1993). 
[9] A. Ghubukov, S. Sachdev, and T. Senthil, J. Phys.: Gondens. Matter 6, 8891 (1994). 
[10] P. Lecheminant, B. Bernu, G. Lhuilher, and L. Pierre, Phys. Rev. B 52, 9162 (1995). 
[11] A. E. Feiguin, G. J. Gazza, A. E. Trumper, and H. A. Geccatto, Phys. Rev. B 52, 15053 
(1995). 

[12] A. E. Trumper, L. O. Manuel, G. J. Gazza, and H. A. Geccatto, Phys. Rev. Lett. 78, 
2216 (1997). 

[13] A. Auerbach and D. Arovas, Phys. Rev. Lett. 61, 617 (1988); Phys. Rev. B 38, 316 

(1988). 

[14] H. J. Schulz, T. A. L. Ziman, and D. Poiblanc, J. Physique I, 6, 675 (1996). 



8 



FIGURES 



FIG. 1: Ground-state energy per site Eq^ for the TLHA as a function of the number of 
sites N in the cluster. Open and full circles give the saddle-point and fluctuation-corrected 
results, respectively; pluses are exact numerical results from and crosses are the rotational 
corrected values from 0|. The line indicates the extrapolation to the thermodynamic limit. 
Inset: Extrapolation to the infinite lattice along the series of clusters with periodic (full 
circles) and twisted (full squares) boundary conditions. 

FIG. 2: Parallel stiffness p\\ for the TLHA as a function of the number of sites in 
the cluster. Open and full circles give the saddle-point and fluctuation-corrected results, 
respectively. The full lines indicate the thermodynamic-limit extrapolations and the dashed 
line is the first-order spin-wave result of |T( 



FIG. 3: The renormalizaton factor Zy = P||/P|j"^ for the TLHA on small clusters. Open 
and full circles correspond to saddle-point (SP) and fluctuation-corrected (FL) results, re- 
spectively. The exact values (EX) and first-order spin-wave (SW) results, indicated by 
crosses and open squares respectively, are taken from . 



FIG. 4: Spin stiffness p\\ for the J1—J2 model as a function of the number of sites A^ in the 
cluster. Top panel: r] = 0; bottom panel: 77 = 0.5. Open and full circles give the saddle-point 
and fluctuation-corrected results, respectively. The lines indicate the thermodynamic-limit 
extrapolation. 

FIG. 5: Spin stiffness py for the J1 — J2 model extrapolated to the thermodynamic limit. 
Dashed and full lines correspond to saddle-point and fluctuation-corrected results, respec- 
tively. Inset: Extrapolation of the critical value r/c = J2c/ Ji where the Neel-order stiffness 
vanishes. 

FIG. 6: Slope a of the scaling Eq^ vs. A^~^/^ on clusters with antiperiodic boundary 
conditions. The inset shows the scaling of the ground-state energy E^^ on clusters with 
periodic (full circles) and antiperiodic (full squares) boundary conditions for rj = 0.5. 
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